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I Abstract 

We discuss the concepts of pseudo-dual frames and approximately 
dual frames, and illuminate their relationship to classical frames. Ap- 
proximately dual frames are easier to construct than the classical dual 
frames, and might be tailored to yield almost perfect reconstruction. 

For approximately dual frames constructed via perturbation the- 
ory, we provide a bound on the deviation from perfect reconstruction. 
An alternative bound is derived for the rich class of Gabor frames, by 
using the Walnut representation of the frame operator to estimate the 
deviation from equality in the duality conditions. 

As illustration of the results, we construct explicit approximate 
duals of Gabor frames generated by the Gaussian; these approximate 
duals yield almost perfect reconstruction. Amazingly, the method 
applies also to certain Gabor frames that are far from being tight. 

1 Introduction 

Let 7i be a separable Hilbert space with inner product (■, ■). Given a frame 
{fk} for H, it is well known that there exists at least one dual frame {gk} , 
i.e., a frame for which 

f = J2^f,fk)gk, \ffen. 

Unfortunately, it is usually complicated to calculate a dual frame explicitly. 
Hence we seek methods for constructing approximate duals. The current 
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paper treats the concepts of pseudo-dual and approximately dual frames and 
examines their properties. The idea of looking at approximately dual frames 
has appeared several times in the literature, for example [U |8l [11] in the 
wavelet setting and [7] for Gabor systems, but our contribution is the first 
systematic treatment. 

An approximately dual frame {gk} associated with {f^} satisfies 

\\f-J2^fJk)9k\\<e\\fl yfeH, (1) 

for some e < 1. (This bound is particularly interesting when e is small.) We 
will show that any approximately dual frame generates a natural dual frame 
in the classical sense; further we obtain a family of frames that interpolates 
between the approximately dual frame and the "true" dual frame. 

We also investigate the use of perturbation theory to construct approx- 
imately dual frames. There are situations where it is hard to find a dual 
frame for a given frame {fk}, whereas a frame {hk} close to {fk} can be 
found for which a dual frame {gk} is known explicitly. We present conditions 
under which such a frame {gk} is approximately dual to {fk}, both in the 
case where {gk} is an arbitrary dual of {hk} and for the particular case where 
it is the canonical dual. 

General Hilbert space estimates might of course be improved in concrete 
cases. For Gabor frames in L^(R) we present a direct way to obtain an 
inequality of the type ([T]), based on the Walnut representation of the frame 
operator [15]; the obtained value e measures the deviation from equality in 
the Gabor duality conditions. 

As an illustration of the method, we construct approximately dual frames 
for two Gabor frames generated by the Gaussian. In these concrete cases, 
we obtain approximations of type ([1]) with very small values for e. 

The paper is organized as follows. Section [2] contains a few elementary 
definitions and results from standard frame theory. In Section [3] we intro- 
duce the new concepts of pseudo-dual and approximately dual frames, and 
investigate their properties. The (apparently) most useful concept of these, 
approximately dual frames, is studied in more detail in Section |H in particu- 
lar, we show how to obtain approximate duals via perturbation theory. The 
direct estimates for Gabor frames are found in Section \5[ and the applica- 
tions to Gabor frames generated by the Gaussian are in Section [61 Finally, 
an appendix collects a few examples and proofs related to classical frame 
theory. 



2 



2 Basic frame theory 



For the introductory frame material that follows, see any standard reference 
on frames, such as [5], [10], [16], [21 Chapter 5]. 

A sequence {fk} in H, indexed by an arbitrary countable index set, is a 
frame if there exist constants A,B > such that 

Any numbers A, B such that ([2]) holds are called (lower and upper) frame 
hounds. 

If at least the upper inequality holds, then {f^} is a Bessel sequence. For 
later use, note that if {fk\ is a Bessel sequence then so is {W fk] for any 
bounded operator W : Ti ^ Ti. 

Given a Bessel sequence {fk}-, the synthesis operator T : C.'^ ^ Ti given 

by 

T{ck} = y^^ckfk 

is linear and bounded, with ||T|| < \fB] the series converges unconditionally 
for all {cfc} G The adjoint of T is the analysis operator T* : Ti. given 
by 

T*f = {{f,fk)}. 

The frame operator is TT*. The frame operator is particularly useful if 
{fk} actually is a frame, but it is well defined and bounded if we just assume 
{fk} is a Bessel sequence. 

Given a Bessel sequence {fk}, a Bessel sequence {gk} is called a dual 
frame if 

f = J2^fJk)9k, ^fen. (3) 

Condition ([3]) means that analysis using {fk} followed by synthesis using {gk} 
yields the identity operator. We adopt the terminology of signal processing 
and speak then about perfect reconstruction. 

One choice of dual frame is the canonical dual frame {{TT*)~^ fk}. When 
{fk} is redundant, infinitely many other dual frames exist. 

Note that if {fk} has the upper frame bound B and {gk} is a dual frame, 
then an application of the Cauchy-Schwarz inequality on 

\\fr = Y.^f,h){gkJ) 
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shows that {g^} satisfies the lower frame inequahty with bound 1/B; thus 
the sequences {fk} and {gk} are indeed frames. For later use, observe that no 
estimate on the upper frame bound for {g^} can be deduced from the duality 
condition or from knowledge about the frame bounds for {fk}, as Example 
lA.ll in Appendix |X] shows. 

3 Various concepts of duality 

In order to apply the dual frame expansion ([3]) for a given frame {fk}, we 
need to find a dual frame {gk}- Unfortunately, it might be cumbersome — 
or even impossible — to calculate a dual frame explicitly. In the literature 
one finds only a few infinite dimensional non-tight frames for which a dual 
has been constructed. This paucity of constructions leads us to seek frames 
that are "close to dual" . We will propose two such concepts. 

Suppose that {fk} and {gk} are Bessel sequences in TC; we will denote 
their synthesis operators by T : £^ 7i and U : i"^ ^ Ti, respectively. The 
operators TU* : H ^ H and UT* : H ^ H will be called mixed frame 
operators. Note that 

Tu*f = J2{f^9k)fk, fen, 

and 

UT*f = J2{fJk)gk, fen. 

Recall that {fk} and {gk} are dual frames when TU* = I or UT* = I. 
Definition 3.1 Bessel sequences {fk} an-d {gk} are said to be 

• approximately dual frames if \\I — TU*\\ < 1 or \\I — UT*\\ < 1, 

• pseudo-dual frames ifTU* or UT* is a bijection on n. 

In each definition, the two given conditions are equivalent, by taking adjoints. 
Note also that the article [1] uses the term approximate dual with a different 
meaning than we do here. 

We illustrate the pseudo-dual idea by means of a frame characterization. 
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Theorem 3.2 (Characterization of frames) Let {f^} be a Bessel sequence 
in Ti. Then the following statements are equivalent. 

(a) {fk} is a frame for Ti. 

(b) The frame operator TT* is a bijection on Ti. 

(c) The synthesis operator T is surjective onto Ti. 

(d) {fk} has a dual frame {gk\- 

(e) {fk} has a pseudo- dual frame {gk}- 

The new part of the theorem is the pseudo-dual statement in part (e). 
We need, however, some information from the proof of the known parts at 
several instances in our paper, so we include a full proof of Theorem 13.21 in 
Appendix El 

Part (b) of the theorem shows that every frame is pseudo-dual to itself. 
The proof of (a)=r-(b) gives even more: in fact, denoting the frame bounds 
for {fk} by A and B, the result f[T6l) in Appendix 1X1 shows that 

2 - - 1 

III-- TT*\\<4 < 1, (4) 

" A + B " - I + 1 ' ^ ' 

i.e., every frame {fk} is approximately dual to the multiple {-^^fk} of itself. 
Next we establish relations among the three duality concepts. 

Lemma 3.3 (Duality relations) Let {fk} one? {gk} be Bessel sequences in 

n. 

(i) If {fk} and {gk} are dual frames, then {fk} and {gk} are approximately 
dual frames. 

(a) If {fk} and {gk} are approximately dual frames, then {fk} and {gk} 
are pseudo-dual frames. 

(Hi) If{fk} and {gk} are pseudo- dual frames, then {fk} and {gk} are frames. 

(iv) If {fk} and {gk} are pseudo-dual frames and W : H ^ H is a bounded 
linear bijection, then {fk} and {Wgk} are pseudo- dual frames. 
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Proof. Statement (i) is an immediate consequence of the definitions. So 
is (ii), since ||/ — UT*\\ < 1 implies UT* is a bijection, with inverse given 
by a Neumann series. Statement (iii) follows from Theorem 13. 2[ For the 
proof of (iv), observe that the synthesis operator for {Wgk} is X = WU; 
the assumptions of {fk} and {gk} being pseudo-dual frames and W being a 
bijection imply that XT* = WUT* is a bijection. □ 

The property of being a pair of pseudo-dual frames is significantly weaker 
than being a pair of dual frames. Nevertheless, the following proposition 
shows that pseudo-dual frame pairs generate dual frame pairs in a natural 
fashion. 

Proposition 3.4 (Pseudo-duals generate duals) If {fk} and {gk} are 
pseudo-dual frames, then {fk} and {{UT*)^^gk} are dual frames. 

Proof. Assuming that {fk} and {gk} are pseudo-dual frames, we know that 
(f/T*)^^ exists and is bounded. Hence {{UT*)^^gk} is a Bessel sequence. If 
we analyze with {fk} and synthesize with {{UT*)^^gk} then we obtain the 
identity, since 

J2{f,fk){UT*)-'9k = {UT*)-'Y.^fJk)gk 
= {UT*)-^UT*f 
= /• 

Thus {fk} and {{UT*)~^gk} are dual frames. □ 

The relation of being a pair of pseudo-dual frames is symmetric: if {fk} 
and {gk} are pseudo-dual frames, then so are {gk} and {fk}- The pseudo-dual 
relation is also refiexive on the set of frames since every frame is pseudo-dual 
to itself, as remarked after Theorem 13.21 Transitivity of the pseudo-dual 
relation fails (in Hilbert spaces of dimension at least 2) by Example 14.41 
below. 

4 Approximately dual frames 

In this section we focus on approximately dual frames {fk} and {gk} . As 
before, we denote the associated synthesis operators by T and U, respectively. 
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Our main goal is to demonstrate how one can use this concept to obtain what 
in engineering terms would be called "almost perfect reconstruction" . 
Recall that if {fk} and {gk} are approximately dual frames then 

11/ - J2^f, fk)9k\\ = \\{I- UT*)f\\ <\\I- UT*\\ \\f\\. 

If ||/ — UT*\\ « 1 then we see that '^{f,fk)gk "almost reconstructs" /, 
which motivates our terminology. 

We first strengthen Proposition I3.4[ for approximately dual frames and 
the associated "natural" dual frame. 

Proposition 4.1 Assume that {fk} o,nd {gk} are approximately dual frames. 
Then the following hold: 

(i) The dual frame {{UT*)~^gk} of {fk} can he written 

oo 

{UT*)''g, = g, + - UTT9k. (5) 

n=l 



(a) Given N ^ N, consider the corresponding partial sum, 

N 

7r = 9k + - UTTgk. (6) 



n=l 



Then {7^^^} is an approximate dual of {fk}- Denoting its associated 
synthesis operator by Z , we have 



I - ZT*\\ < \\I - UT*\f^\ (7) 



Proof. Assuming that {fk} and {gk} are approximately dual frames, the 
inverse of UT* can be written via a Neumann series as 

00 

{UT*y^ =(/-(/- UT*)) = J](J - f/T*)". (8) 

n=0 

The result in (i) now follows by applying this expansion to gk and recalling 
Proposition 13.41 
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For (ii), note that {7^^''} is a Bessel sequence since it is obtained from 
the Bessel sequence {gk} by a bounded transformation. And 

N 

ZT*f = J2{fJk){l + J2{I-UTT)9k 

n=l 

N 



- UT*YUT*f 



n=0 

N 



= UT*Y [I -{I- UT*)) f 

n=0 

= f - {I -UTY^^f 

by telescoping. Thus 

\\I-ZT*\\ = \\(I -UT*)^+^\\ 



< \\I -UT*\f+^ <1. 

□ 



Remark 4.2 In view of formula in Proposition 14. 1^ an approximately 
dual frame {gk} associated with a frame {fk} can be regarded as a zero- 
th order approximation to the (exact) dual frame {{UT*)^^gk}. In case 
||/ — UT*\\ is small, reconstruction using the approximate dual {gk} is close 
to perfect reconstruction. The result in (ii) yields a family of approximately 
dual frames that interpolates between the approximate dual {gk} and the 
dual frame {{UT*)~^gk}; the estimate ([7]) shows that by choosing suf- 
ficiently large, we can obtain a reconstruction that is arbitrarily close to 
perfect. The drawback of the result with respect to potential applications is 
the complicated structure of the operator in ([6]) defining the sequence 7^'^''. 

We have now presented the basic facts for approximately dual frames. In 
the following, we will present concrete settings where they yield interesting 
new insights. 

To motivate the results, we note (again) that it can be a nontrivial task 
to find the canonical dual frame (or any other dual) associated with a general 
frame {fk}- We seek to connect this fact with perturbation theory by asking 
the following question: if we can find a frame {hk} that is close to {fk} and 
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for which it is possible to find a dual frame {gk}, does it follow that {gk} is 
an approximate dual of {/fc}? We will present some sufficient conditions for 
an affirmative answer. 

First we state a general result, valid for dual frames {gk} with sufficiently 
small Bessel bound; later we state a more explicit consequence for the case 
where {gk} is the canonical dual frame, see Theorem l4.5[ Technically, we need 
not assume that {fk} is a frame in the general result: the frame property 
follows as a conclusion. On the other hand, as explained above, the main use 
of the result is in a setting where {fk} is known to be a frame in advance. 

Theorem 4.3 (Dual of perturbed sequence) Assume that {fk} is a se- 
quence in Ti and that {hk} is a frame for which 

Y.\{fJk-hk)\'<R\\f\\\ yfen, 

for some R> 0. Consider a dual frame {gk} of {hk} with synthesis operator 
U, and assume {gk} has upper frame hound C . 

IfCR < 1 then {fk} o-nd {gk} are approximately dual frames, with 

\\I - UT*\\ < \\U\\ < \fCR < 1. 

Proof. With our usual notation, we have UV* = I since {gk} and {hk} are 
dual frames. Hence 

||/ _ UT*\\ = \\U{V* - T*)\\ < \\U\\ \\V* - T*\\ < VCR < 1. 

□ 

It is crucial in Theorem 14.31 that the dual frame {gk} has upper frame 
bound less than 1/R. Otherwise {gk} need not be pseudo-dual to {fk}, let 
alone approximately dual, as the next example shows. 

Example 4.4 Consider = with the standard basis {61,62}. Let e > 
and consider the frames 

{fk} = {0, 61, 62}, {hk} = {eei, 61, 62}, {gk} = {e"^6i, 0, 62}. 

Write T, V, U for the associated synthesis operators. Note that 

TU*f={f, 62)62; 
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this operator is neither injective nor surjective, and so {fk} and {g^} are not 
pseudo-dual frames. Clearly, {h/.} is a frame for C^, regardless of the choice 
of e. Now, because 

3 
k=l 

the condition in Theorem 14.31 is satisfied with R = e^. Thus we see that no 
matter how close {hk} gets to {fk} (meaning, no matter how small e is), 
the frame {gk} is not pseudo-dual to {fk}- Note that for e < 1, the frame 
{gk} has the upper frame bound C = = 1/R; thus. Theorem 14.31 is not 
contradicted. 

In this example it is immediate to see that TV* = I and VU* = /; thus 
{fk} and {hk} are dual frames, and so are {hk} and {gk}- In particular, the 
example shows that the pseudo-dual property is not transitive. The same 
argument applies in any separable Hilbert space of dimension at least 2. □ 



For practical applications of Theorem 14. 3 [ it is problematic that the Bessel 
bound C on {gk} might increase when one reduces R by taking {hk} very 
close to {fk} : in this way, it might not be possible to get CR < 1 just by 
taking {hk} sufficiently close to {fk}- In the case when {gk} is the canonical 
dual of {hk} the problem does not arise: as the next result shows, if {hk} 
is sufficiently close to {fk}, then the canonical dual {gk} of {hk} is approx- 
imately dual to {fk}- For the statement of this result we need to assume 
that {fk} is a frame, since its lower frame bound controls the perturbation 
hypothesis. 

Theorem 4.5 (Canonical dual of perturbed frame) Let {fk} be a frame 
for Ti with frame hounds A, B. Let {hk} he a sequence in Ti for which 

Y.\{fJk-hk)\'<R\\f\\\ yfen, 

for some R < A/ A. Denote the synthesis operator for {hk} hy V. 

Then {hk} is a frame. Its canonical dual frame {gk} = {{VV*)~^hk} is 
an approximate dual of {fk} with 

11/ - UT*\\ < , L < 1, 

where U denotes the synthesis operator for {gk}- 
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Proof. The sequence {h^} is a frame with frame bounds {^/A — ^/R)"^ and 
(V^ + VR^, by [2, Corollary 5.6.3]; the proof is short, and so we give it 
here. We have ||T* — V^*||^ < -R by assumption, and so 

irvii^^- 11^711^^1 <v^ii/ii 

for all / G 7i, from which the desired frame bounds follow easily. This part 
of the proof uses only R < A. 

The canonical dual frame of {hk} is {gk} = {{VV*)~^hk}, with frame 
bounds 

and 



by [21 Lemma 5.1.6]. In terms of the synthesis operator U for {gk}, the upper 
bound says 

1 

< 



y/A-y/R' 
Theorem 14.31 therefore implies 

\\I-UT*\\ < \\U\\y/R 



< 



R 



^fA-\fR 
1 

To complete the proof, just notice this last expression is smaller than 1 if 
and only \iR< Aj^. □ 

We end this section with remarks concerning the appearance of approxi- 
mate duals in the literature. Approximately dual frames have been employed 
in the theory of wavelets; for example, Holschneider developed a suffi- 
cient condition for a pair of wavelet systems to be approximately dual [HI 
Section 2.13], and Gilbert et al. [S] obtained approximate duals for highly 
oversampled wavelet systems by perturbing the continuous parameter sys- 
tem. Approximate duals have proved central to the recent solution of the 
Mexican hat wavelet spanning problem [1]. There the synthesizing wavelet 
frame is given, being generated by the Mexican hat function, and the task 
is to construct an approximately dual analyzing frame. Interestingly, this 
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approximate dual construction holds in all L'^ spaces (in the frequency do- 
main), which suggests that parts of the theory in this paper might extend 
usefully to Banach spaces. 

Approximately dual frames in Gabor theory arose in the work of Fe- 
ichtinger and Kaiblinger [71 Sections 3,4]. They studied the stability of Gabor 
frames with respect to perturbation of the generators and the time-frequency 
lattice; among many other results, they showed that a sufficiently small per- 
turbation (measured in a specific norm) of a dual frame associated with a 
Gabor frame leads to an approximately dual frame. In the next sections, we 
construct explicit examples of approximately dual Gabor frames. 

5 Gabor frames and approximate duals 

The general frame estimates in Section H] are based on operator inequalities 
such as the product rule \\TU\\ < \\T\\\\U\\. We will apply these general 
estimates in the next section to obtain explicit approximate duals for Gabor 
frames generated by a Gaussian. First, though, we will sharpen the general 
frame estimates in the concrete, rich setting of Gabor theory. These sharper 
estimates also will be applied in the next section. 
A Gabor frame is a frame for L^(R) of the form 

(e2--^^(x-na)}„,„gz 

for suitably chosen parameters a,b > and a fixed function g G L^(M), called 
the window function. The number a is called the translation parameter and 
b is the modulation parameter. Introducing the operators 

{Tag){x) = g{x - a), {E,g){x) = e^^''"'g{x), for a, 6, x G M, 

the Gabor system can be written in the short form {EmhTna9}m,n&- For 
more information on Gabor analysis and its role in time-frequency analysis 
we refer to the book by Grochenig [9]. 

We state now the duality conditions for a pair of Gabor systems, due to 
Ron and Shen We will apply the version presented by Janssen |12j : 

Lemma 5.1 Two Bessel sequences {EmbTn^}m,n& and {EmbTng}m,n& form 
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dual frames for if and only if the equations 

''^^(p{x — ak)g{x — ak) — b = 0, (9) 

A;GZ 

ip{x - n/b - ak)g{x - ak) = 0, nGZ\{0}, (10) 



fcez 



hold a.e. 



Recall that the Wiener space W consists of all bounded measurable func- 
tions : M — i> C for which 



Y\\gX[n,n+l[\\oo < OO. 



It is well known that if g E W then {EmbTnag}m,nei is a Bessel sequence for 
each choice of a, 6 > 0. 

Theorem 5.2 Given two functions ip,g E W and two parameters a,b > 
0, let T denote the synthesis operator associated with the Gabor system 
{EmbTna^}m,neii ^^^d U the synthcsis operator associated with {E.mbTnag}m,nei 
Then 



\I -UT*\\ < 



1 



b-'Y] TakgTakV + Tn/bTakgTakV 



k£Z 



Proof. The starting point is the Walnut representation of the mixed frame 
operator associated with {EmbTna^]m,n& and {EmbTnag]m,n(ii- According to 
Theorem 6.3.2 in [9J, 

UT*f{-) = \Y. [Y.T,kV{--n/b)T,,g{-)\ T^/bf{-). f G L^R). 
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Thus, 



\\f~UT*f\\ 

- i:h -^'^o^kfTakg WfW + ^Y\\YTn/bTakfTakg 



k£Z 



nj^O k& 



which concludes the proof. 



□ 



Observe that the terms appearing in the estimate in Theorem 15.21 mea- 
sure the deviation from equahty in the duahty relations in Lemma 15. with 
respect to the ||-||oo-norm. In particular. Theorem 15.21 proves the "sufficient" 
direction of Lemma 15.11 under the stated hypotheses. 



6 Applications to Gabor frames generated by 
the Gaussian 

The Gaussian g{x) = e~^^ is well known to generate a Gabor frame when- 
ever ab < 1. For the case where the parameters a and b are equal and small, 
Daubechies has demonstrated that the Gabor frame generated by the Gaus- 
sian is almost tight, see [Bj, p. 84-86 and [6j, p. 980-982; in particular, the 
formula (jlj) yields almost perfect reconstruction. 

Regardless of the choice of a and b, no convenient explicitly given ex- 
pression is known for any of the dual frames associated with the Gaussian, 
though. We will construct explicit, approximately dual frames associated 
with the Gaussian, for certain choices of a and b. These approximately dual 
frames provide almost perfect reconstruction; Example 16. 21 is particularly in- 
teresting because it deals with a frame that is far from being tight, i.e., no 
easy way of obtaining an approximately dual frame is available. 

Example 6.1 Consider the (scaled) Gaussian 

151 



'111 
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The scaling is introduced for convenience; similar constructions can be per- 
formed for other Gaussians as well. It is well known that for this Gaussian, 
the Gabor system {EmbTnip}m,n& forms a frame for L^(R) for any suffi- 
ciently small value of the modulation parameter 6 > 0; we fix 6 = 0.06 in 
this example. Denote the synthesis operator for that frame by T. 

We will use the results derived in this paper to find an approximately dual 
Gabor frame. Let h = denote the eighth-order B-spline, centered at the 
origin. This B-spline approximates the Gaussian very well: see Figures [1] 
and [21 The Gabor system {Ei,rnTnh}m,n£i is a frame for L^(M), and by 
Corollary 3.2 in [3] it has the (non-canonical) dual frame {EbmTng}m.,n£Z, 
where 

7 

g{x)=bJ2 M^ + n) (12) 

n=-7 

is a linear combination of translated B-splines. This g is supported on 
[—11, 11] and is constantly equal to b on supping = [—4,4]; see Figure [31 

We first apply Theorem 14.31 A numerical calculation based on Theorem 
9.1.5 in [2j| shows that the functions {EbrnTn{(p — h)}m,n&z form a Bessel 
sequence with Bessel bound R = 0.0006. 

Furthermore, a second application of Theorem 9.1.5 in [2j reveals that 
{EhmTng}m,n€Z has Besscl bound C = 1. Thus the hypotheses of Theorem 
14.31 are satisfied. Denoting the synthesis operator for {-Ef,m^n5'}m,nez by U, 
we deduce that 

\\I-UT*\\ < VCR< 0.025. 

On the other hand, an application of Theorem 15.21 immediately yields the 
better estimate 

\\I -UT*\\ < 0.0027. 
In terms of the Gabor systems involved, this conclusion says that 
II/- J2 {f,E,mTnV)EbmTng\\ < 0.0027 || / || 2, V/ G L\R). 

That is, analysis with the scaled Gaussian and synthesis with the function g 
yields almost perfect reconstruction. 
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Figure 3: The non-canonical dual window g of Bg, see (|T2l) . 

We note that the frame considered in this example is almost tight; thus, 
the approximate dual constructed here does not perform better than using a 
scaled version of the frame itself as approximate dual. For a frame that is far 
from being tight this aspect changes drastically in favor of the approximate 
duals considered in this paper - see Example I6.2[ Lastly, we remark that 
other approximations to the Gaussian are possible too; instead of using B- 
splines and their duals from ^ we could have used certain other splines and 
their duals from [T3] . 

□ 

In the next example we apply Theorem 14. 51 and Theorem l5.2[ Let us again 
consider B-spline approximation to a Gaussian. The advantage of using the 
canonical dual and Theorem 14.51 compared with our use of a non-canonical 
dual frame in Example 16. ![ is that larger values of the modulation parameter 
b can be handled. Approximating using Bg, our approach in Example 16.11 is 
restricted to b < 1/15 (the restriction comes from the underlying results in 
[3]); as illustration of the larger range for b obtained via Theorem 14.51 we 
take 6 = 0.1 in the next example. 

Example 6.2 Let 
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An application of Theorem 5.1.5 in [2] shows that the functions {-E'o.im^nV^jm.nGZ 
form a frame for L^(M) with frame bounds A = 2.6, B = 10.1. The function 

315 

h{x) = —Bs{2Mx) 
151 

yields a close approximation of (p, see Figure 4. A numerical calculation based 
on Theorem 5.1.5 in [2] shows that the functions {Eo imTn^f — h)}m,nGZ form 
a Bessel sequence with Bessel bound R = 6.5 ■ 10"^ < A/A. 

The function h has support on [—4/2.36,4/2.36], an interval of length 
8/2.36. Since the modulation parameter b = 0.1 is smaller than (8/2.36)"^ 
and the function 

H{x) = J2\Hx + k)\'' 

is bounded above and below away from zero, it follows from Corollary 9.1.7 in 
[2j that the frame {Eo.imTnh}„,^n& has the canonical dual frame {-E'o.im^nfi'jm.nGZ, 
where 



15.1 1 



Bs{2.36x). (13) 



315 E„GJ^8(2.36(x + n))P 
See Figure 5. 

The frame {Eo,i„^T„g}m,nez is approximately dual to {Eo,imTn(p}m,n&, by 
Theorem 14.51 Denoting their synthesis operators by U and T respectively, 
the approximation rate in Theorem 14.51 is measured by 

11/ _ < — i < 0.016. (14) 

- ^/A/R-1 - 

On the other hand. Theorem 15.21 yields the somewhat better estimate 

\\I-UT*\\ < , ^ < 0.009. (15) 

Thus the approximate dual frame almost yields perfect reconstruction. 

Note that the frame {-Eo.im^nV^}m,nGZ is far from being tight. Thus, the 
sequence {i^o.im^n (34^ V^) }m,nez is a poor approximate dual: the estimate 
in (jlj) yields 

9 B _ 1 

S\\ < ^—^ = 0.59, 



A + B " - f + 1 
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Figure 5: The approximative dual window g in ( |T3i) . 
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which is far worse than the result in (fT^ . 



□ 



Remark 6.3 

1. In closer approximation than in (fT5|) is required, an application of 
the iterative procedure in Proposition l4.1l (ii) can bring us as close to perfect 
reconstruction as desired. Letting {fk} be a Gabor frame {EmhTnaf}m,n& 
and {gk} be an approximately dual Gabor frame {EmbTnag}m,n&-i each ap- 
proximate dual {7^^^} constructed in Proposition 14. II will again have Gabor 
structure. For = 1, an easy calculation shows that {7^^"*} will be the 
Gabor system {E„,bTnal}m,n&, where 

-i = g+{I-UT*)g = 2g-UT*g 

= 2g — {g, Em'bTn'af)Em'bTn'ag 

m',n'£Z 

= C^- {9J))9- ^ {g,E^>i,Tn'af)E„^'bTn'a9- 
(m',n')7^{0,0) 

For the particular case considered in Example 16. 2[ and denoting the synthesis 
operator for {EmbTnal}m,n&z by Z, Proposition 14.11 shows that the estimate 
in f fT4|) will be replaced by 

||/ - Z*T\\ < 0.009^ = 8.1 X 10~^ 

2. Another way of calculating approximate duals for Gabor frames based 
on (p{x) = e~^^ would be to follow the procedure in Example 16.21 with a 
function h of the type 

h{x) = (p{x)xiix) 

for a sufficiently large interval /, rather than letting hhe a B-spline. However, 
for such a function h the canonical dual window associated with {-Eo.im^nV^}m,nGZ 
will not be continuous; thus, it will have bad time-frequency properties. 
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A Proofs and additional examples 



In this Appendix we collect some additional information related to classical 
frame theory. 

First, we present the example announced in Section [2], showing that an 
estimate on the upper frame bound for a frame {gk} neither can be deduced 
from the fact that it is dual to a frame {gk}, nor from knowledge of the frame 
bounds for {fk}- 

Example A.l Consider the Hilbert space 7i = with the standard or- 
thonormal basis {61,62}. Let 

{/l,/2,/3} = {0,61,62}. 

Then {fk} is a frame with bounds A = B = 1. For any C G M, the sequence 

{5-1, 5-2, 5-3} = {^61,61,62} 

is a dual frame of {fk}- The (optimal) upper frame bound for {^f^} is + 1, 
which can be arbitrarily large. 

Exactly the same considerations apply in an arbitrary separable Hilbert 
space of dimension at least two. □ 



We now state the announced proof of Theorem 13. 2[ The equivalences 
between the first four statements are well known; we include the proof in 
order to keep the paper self-contained, and because we refer to some of the 
steps elsewhere in the paper. 

Proof of Theorem 13. 2t First recall that T and U are bounded operators 
from £^ to H, whenever {fk} and {gk} are Bessel sequences. 

(a) implies (b): First, write S = TT* for the frame operator, and note it 
is selfadjoint. By part (a) we have < (5/, /) < B\\ff for all / G H. 

That is, AI < S < BI; following the argument by Grochenig [9] p. 91, we 
inter that -S^I < I - < ^I. Thus, 

< 1. 

Therefore -j^S is invertible, so that S* is a bijection as desired. 
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(b) implies (c) is immediate. 

(c) implies (d): consider the pseudo-inverse -.n f. Applying the 
Riesz representation theorem to the functional {T^^ f)k = (/c-th component 
of T'f/), we obtain some element gk & H satisfying {T^f)k = {f,gk)- The 
sequence {gk} is Bessel, since 

J2\(f^9k)\' = Yl = iitvii,'. < iiTi,^.ii/f . 

The analysis operator U* associated to the sequence {g^} equals exactly the 
pseudo-inverse T\ by construction, and so TU* = TT^ equals the identity, by 
the fundamental property of the pseudo- inverse. Thus {gk} is a dual frame 
for {/fc}, and so (d) holds. 

(d) implies (e) trivially. 

(e) imphes (a) as follows: suppose {gk} is a pseudo-dual frame for {/fc}, 
so that TU* is a bijection on Ti. Then for all f ^Ti, 

\\f\\ = \\{TUT'TU*f\\ 

< \\{TU*)-'T\\\\U*f\\i^, 

and of course ||f/*/||^2 < ||f/*|| ||/||. Hence {fk} is a frame with bounds 
A = 1/ \\{TU*)-'^Tf and B =\\U*f. □ 
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